Cold atoms in optical lattices: a Hamiltonian ratchet 
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We investigate a new type of quantum ratchet which may be realised by cold atoms in a double- 
well optical lattice which is pulsed with unequal periods. The classical dynamics is chaotic and we 
find the classical diffusion rate D is asymmetric in momentum up to a finite time t r . The quantum 
behaviour produces a corresponding asymmetry in the momentum distribution which is 'frozen- 
in' by Dynamical Localization provided the break-time t* > t r . We conclude that the cold atom 
ratchets require Db/h ~ 1 where 6 is a small deviation from period-one pulses. 
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Cold atoms in optical lattices provide an excellent ex- 
perimental demonstration of the phenomenon of Dynam- 
ical Localization jlj. Dynamical Localization (DL) has 
been described as the so-called 'quantum suppression of 
classical chaos'. More precisely, in the usual realizations, 
a periodically-driven or kicked system is associated with 
classically chaotic dynamics for sufficiently strong pertur- 
bation. Its average energy grows diffusively and without 
limit, whereas for the corresponding quantum system the 
diffusion is suppressed after an /i-dependent timescale, 
the 'break-time' t*. The final quantum momentum dis- 
tribution is localized with a characteristic exponential 
profile. The formal analogy established with Anderson 
Localization 0] forms a key analysis of this phenomenon. 

A series of recent experiments on Cesium atoms || 
gave a classic demonstration of this effect for atoms in 
driven or pulsed lattices of sinusoidal form. The possi- 
bility of experiments with asymmetric lattices, in par- 
ticular with asymmetric double wells || |5|, leads us 
to investigate the possibility of constructing a 'clean' 
atomic ratchet, where the transport results purely from 
the chaotic Hamiltonian dynamics, with no Brownian or 
dissipative ingredients. 

There is already an extensive body of work on Brow- 
nian ratchets [||, driven by the need to understand bio- 
physical systems such as molecular motors and certain 
mesoscopic systems. Some of this work encompasses the 
quantum dynamics R|. However, to date there has been 
very little work on Hamiltonian ratchets. One notable 
exception is the work by Flach et al || where the gen- 
eral form of the spatial and temporal desymmetrization 
required to generate transport was investigated. To our 
knowledge, the only substantial study of quantum Hamil- 
tonian ratchets to date, however, is the work of Dittrich 
et al || which concluded that transport occurs in mixed 
phase-spaces. They demonstrated that overall transport 
(in their system with periodic boundary conditions) is 
zero if starting conditions cover all regions of phase-space 
(including stable islands and chaotic regions) uniformly. 
A key result was a sum rule which exploited the fact 
that transport in the chaotic manifold is balanced by 
transport in the adjoining regular manifolds (stable is- 
lands/tori). Very recenty 10 , it was shown that a kicked 
map with a 'rocking' linear potential leads to confinement 



in the chaotic region between a pair of tori which are not 
symmetrically located about p = 0. 

Here we propose a new type of Hamiltonian quantum 
ratchet which, classically, is completely chaotic. This 
ratchet is consistent with the rules established by Dittrich 
et al , but relies on a different mechanism: the transport 
is due to a non-uniform momentum distribution within a 
completely chaotic phase-space. To our knowledge, this 
is the only example of a clean, non-dissipative ratchet 
which does not rely on specific dynamical features such 
as a particular set of islands/tori. 

The basic mechanism is as follows: for a repeating cy- 
cle of kicks, of strength K e ff, perturbed from period-one 
by a small parameter 6, we find that the classical diffu- 
sion rates for positive and negative momenta (D + and 
D~ respectively) are different up to a finite time. Up 
to this time, an asymmetry in the classical momentum 
distribution N c i(p) accumulates with kick number. Be- 
yond this 'ratchet' time, t r , the rates equalize, we have 
D + ~ D~ ~ D (where D ~ K* f J is the total diffusion 
rate) , no more asymmetry accumulates and the net classi- 
cal momentum < p c i > becomes constant. But < p 2 cl >, 
of course, continues to grow with time as ~ Dt. We 
find that the corresponding quantum current depends on 
t*/t r : if the break time is too short no asymmetry in the 
quantum N qm (p) accumulates and there is no quantum 
transport. If t* >> t r the localization length L becomes 
large and the effective quantum momentum asymmetry 
~< Pqm > /L decreases. We find that t r ~ A 
quantum ratchet will have the clearest experimental sig- 
nature if t* ~ t r . Since t* ~ D/h 2 , our main conclusion 
is that optimal cold atom ratchets need Db/h ~ 1 with 
small b < 0.2. 

Here we consider a system with a Hamiltonian given by 
H = ^ + KV{x) £ n> . S(t - nTi). We model the double- 
well ratchet potential as V(x) = {sinx+asin(2a;+$)}. In 
the usual realization of DL, the Quantum Kicked Rotor 
(QKR), the pulses or kicks are equally spaced. For the 
ratchet, we introduce a repeating cycle of N unequally 
spaced kicks, ie T\, T 2 ..Tn. 

The time-evolution operator for the i — th kick of the 
n — th cycle factorises into a free and a 'kick' part Ui — 

jjfreejjkick 
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FIG. 1: Effect of equal and unequal kick spacings on a mini- 
mal uncertainty gaussian wavepacket with h — 1/2 with ini- 
tial < p >— 0, for kick strength K = 2 and a = 1/2, for 
different 3>. (a) Shows evolution of energy and momentum 
for T=l (dw-QKR) (b)Shows evolution of energy and mo- 
mentum for a repeating AT = 3 cycle of kicks with b — 0.1, 
hence Ti = 1.1,1.0,0.9 (cdw-QKR). Figs (a) and (b) show 
that the quantum break time is substantially longer for the 
unequal kicks. In (b) we illustrate our finding that the classi- 
cal diffusion rate in the chirped case is very well fitted by 
the lowest order approximation Do ~ K%ff/2 ~ 4 where 
K 2 ff = K 2 (l + 4a 2 ) for D ~ 2 or larger. They also show 
that there is no transport for the equal kick case, but that 
there is a substantial net momentum (~ est if t > t*) for the 
cdw-QKR. Figs, (c) and (d) show corresponding quantum 
momentum distributions for the dw-QKR and the cdw-QKR. 
The figures show that while the logarithmic distribution takes 
the usual triangular form for the dw-QKR, for the cdw-QKR 
an asymmetric distribution results. In the lower figure we 
plot the first moment of the momentum distribution: (to ease 
comparison we use \p\). The net momentum is the difference 
in area between the positive p and negative p regions. The DL 
form \p\.N(\p\) = ^exp-|p|/L (with L = 27.5 ~ S.5D /h 
(see fig. 3)) is superposed, illustrating that the DL form gives 
excellent agreement for all cases for large enough p. 



In the usual plane wave basis, for a given quasi- 
momentum q, the matrix elements of Ui can be shown 
to take the form: 



(0.1) 



where the J are ordinary Bessel functions. The time 
evolution operator for one period U(T = J^. T) = 



is an effective H e ff = 8w r T where uj r is the recoil fre- 
quency. In ||, h e ff ~ 1, so here we have considered 
the range H = 1 — > 1/10. We can use a re-scaled time 
such that, without loss of generality, we take the aver- 
age over the periods in each cycle to unity < T >= 1. 
We define our 'chirped' kicks to mean that we have 
a repeating cycle of kicks separated by time intervals 
1 + rib, 1 + (n — 1)6..., 1 — (n — 1)6, 1 — nb where n > is 
an integer and 6 is a small time increment. N = 2n + 1 
for N odd and N — 2n for N even. In particular, our 
N = 3 cycle corresponds to a repeating set of kick spac- 
ings Ti = 1 + 6, T 2 = 1, T 3 = 1 - 6 while an N = 2 cycle 
corresponds to T\ = 1 + b, T-x = 1 — b and so forth. 

In Fig.l we compared the evolution of a quantum 
wavepacket with equal kick times (Tj = 1) with a cor- 
responding unequal-kick case with N = 3, 6 = 0.1. Since 
V{x) in general represents a double well potential, we 
refer to it as the dw-QKR to distinguish it from the 
standard map case with V = Ksinx. For convenience 
and by analogy with optical terminology, we refer to the 
unequal-kick case as the 'chirped' or cdw-QKR. N = 3 is 
the minimum number needed to break time-reversal in- 
variance: N = 2 in this system gives no transport. Larger 
numbers (JV > 3) do give transport but have the effect 
of increasing break-times, which, as discussed below, we 
wish to keep small. 

The upper graph in Fig.l shows that, in both cases, the 
average energy of a classical ensemble of particles grows 
linearly with time t , ie < p 2 >~ Dt. Neglecting all cor- 
relations, it is approximated by D ~ Dq 



Kl f j/2 where 



Ui. In the experiments, an important parameter 



K e ff — KyJ (1 + 4a 2 ). The <!> dependence appears in ne- 
glected correlations which in this case appear as products 
of Bessel functions [[fl] . For the standard QKR the cor- 
rections appear as oscillations in D(K), with the maxima 
corresponding to the K's for which accelerator modes are 
important . They were observed experimentally in p2| . 
We find that these corrections are unimportant for the 
cdw-QKR : the lowest order approximation gives a good 
fit from quite low values of K. The equal-kick case, the 
dw-QKR does need substantial corrections. In fact, since 
for the cdw-QKR phase-space periodicity in momentum 
is largely destroyed and indeed we find no evidence of 
accelerator modes; the dw-QKR is periodic in p and we 
find strong effects due to accelerator modes. 

The quantum energy follows the classical trend but 
then saturates to a constant energy when < p 2 >^ L 2 , 
where L is the localization length, after the break-time t* . 
The break-time is substantially longer however (typically 
~ 3 times) in the chirped case relative to the dw-QKR. 
We find that t* grows as N. The figure shows that there 
is no net momentum in the T — 1 case and, asymptot- 
ically, < p >~ 0. However, for the chirped case, for 
t > t*, in general we have < p >^ est. A meaningful 
way to quantify the asymmetry is a re-scaled momen- 
tum pl =< p > / L which also tends to a constant for 
t > t* (eg p L ~ 1/8 in Fig. lb for $ = 0.) Taking 
$ = 7r reverses the symmetry of V(x) and the direction 
of motion relative to $ = 0. Intermediate values of $ 
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typically give < p > within these extremes. The exper- 
iments in H measured the momentum distribution of 
a cloud of atoms and observed the hallmark of DL : a 
triangular form of the logarithm of the momentum dis- 
tribution N(p). Fig 1(c) (upper) shows this is also the 
case for the dw-QKR. However, for the chirped case, the 
corresponding distribution is evidently asymmetric. 
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FIG. 2: Behaviour of the cdw-QKR. The results were ob- 
tained from quantum wavepackets and classical trajectories. 
Both the classical and quantum < p > reach a constant value 
after a finite time (t* and t r respectively), (a) Graph show- 
ing that the quantum ratchet transport is very sensitive to 
ft if the quantum wavepacket localizes before the classical 
wavepacket reaches t r . The values of ft, (0.8 — > 0.125) are 
indicated. K = 1.6, a = 1/2, N = 3 (hence D ~ 2.5) and 
b — 0.1. The asymptotic quantum < p > increases with ft 
and 'catches' up with the classical results for ft ~ 0.25. (b) 
Evolution of < p > for a classical wavepacket (500,000 par- 
ticles with a gaussian random distribution in x,p, of width 
a = 1.5 for b = 0.05 but different D. The classical net mo- 
mentum increases and then saturates after a classical satu- 
ration time t r ~ l/(Db 2 ). (c) The effect of different classi- 
cal diffusion rates (D + (t,b), D~ (t,b)): < p 2 > is evaluated 
separately for positive and negative momenta for a cloud of 
classical particles with D ~ 2.5 and different b. We see that 
the growth in < p 2 > shows a divergence from linear growth 
which, for small by < p 2 >, is similar in magnitude but oppo- 
site in sign for the negative and positive components. The + 
and — indicate < p 2 > + , < p 2 >~ respectively. One sees that 
if t > t r then D + ~ D~ ~ D ~ 2.5. The lower graphs show 
< p 2 > (_) - < p 2 > (+) (for these parameters, D~ > D + ) 
and show the t r oc 1/6 2 dependence. 

The form for the first moment of Af(|p|) shows most 
clearly the source of the non-zero average of p. The dw- 



QKR is symmetric. The cdw-QKR distribution, shows a 
substantial asymmetry. These asymptotic 'frozen' distri- 
butions are insensitive to changes in the starting position 
in x. We have investigated this effect for a range of dif- 
ferent K 1 h and a. 

In Figs. 2a and 2b we show both quantum and clas- 
sical net momenta < p > increase in magnitude, then 
saturate to a constant value after a finite time ~ t* in 
the quantum case or a 'classical ratchet' timescale t r in 
the classical case. But they show other very interest- 
ing and surprising features. The quantum dependence 
on ft is shown in Fig. 2a. The < p qm > are negligible 
for h > 1 but increase rapidly with decreasing ft, up to 
ft ~ 0.25. This is important for any experiment : for 
these parameters, (D ~ 2.5, b = 0.1) an experiment with 
h e ff ~ 0.8 would show little asymmetry, but just halv- 
ing H e ff to ~ 0.4 would show substantial asymmetry. 
Beyond h ~ 0.4, < p qm > is comparable to the satu- 
rated classical value. But since the most experimentally 
'detectable' ratchet is one which maximizes the asymme- 
try in Fig. 1(d), this means maximizing a re-scaled mo- 
mentum pl —< Pqm > / L, so there is no advantage in 
reducing h much below ~ 0.4 since L ~ h . 

Fig 2b shows that, for a given b (b = 0.05 in this graph), 
the classical saturation time falls with increasing D. The 
fact that < pd > saturates at all is surprising: after 
all, the ensemble of classical trajectories is continually 
expanding and exploring new phase-space regions corre- 
sponding to higher momenta. Hence specific dynamical 
features such as partial barriers (cantori) do not account 
for this ratchet effect, though we find 'scars' of cantori 
do affect the detail of N(p). 

We solve the classical dynamics of the dw-QKR with 
the usual map , Xi+i = Xi + Pi ; = Pi — V'(xi+i) for 
the i — th kick. For the chirped case, however, we have 
a 3-kick map ( we consider N = 3 here but the analysis 
for N > 3 is qualitatively similar). We consider the first 
kick of the n — th cycle: 

x nl =x n0 + P n0 (l + b) (0.2) 



Pnl = PnO - V'{ Xn0 + P„ (l + b)) (0.3) 

The effect of the chirping is to allow the free evolution to 
proceed for an additional distance S± = P n ob. The second 
kick has corresponding 5% — and the third kick has 
5.3 = —Priib- Since b is a small number, this represents 
a small perturbation at the start of the evolution of the 
cloud and to first order, for the first kick, 

Pnl = PnO - {V'{x n0 + P n0 ) + Pn0bV"(x n0 + P nQ )) 

(0.4) 

In general, if we work out the change in < p 2 > for 
successive kicks in the standard map we obtain a diffu- 
sion rate D which is the same (after the first few kicks) 
whether we average the momentum from — > oo or from 
— > — oo. Both the zero-th order contribution K 2 /2 as 



well as higher corrections (such as those, for example, 
involving correlations between the i — th kick and the 
(i + 2)th kick) are even in p and insensitive to the sign 
of p. However, in the cdw-QKR case, if we can make the 
perturbative expansion above, we now have correlations 
which depend on the sign of p and which scale with b. 
For our classical ratchets, we calculated -separately- 

< p 2 >'~~) for those particles with p < and < p 2 >( + ) 
for those with p > 0. The results are shown in Fig. 2c 
for Do ~ 2.5 and different b. They are quite striking : 

< p 2 >( - ' and < p 2 >( + ' separate gradually, more or 
less symmetrically, about the line ~ 2.5i, but beyond a 
certain time they run parallel to each other and their 
slopes become equal with D^ + ' ~ ~ 2.5. In the 
lower part of Fig. 2c we plot the actual difference < 
p 2 >( _ ) — < p 2 >W for each b since this shows the 
saturation effect more clearly. 

The saturation time t r is very important since then 
the classical ratchet speed reaches its maximum and 
for t > t r , < Pd >^ est. We identify it as a point 
far enough from where we can make the small angle 
perturbative approximation above. This will happen 
when by/ < p 2 >(+) ~ n for the positive component and 
by < p 2 >(~) ~ 7r for the negative component. For an 
order of magnitude estimate of the mean timescale in- 
volved, we take bVDt ~ ir. Hence we obtain t r ~ 

Numerically, we estimate t r ~ which is not in- 
consistent with the above. This explains the counter- 
intuitive behaviour that the larger deviation from period- 
one kicking (ie the larger b), give a smaller ratchet effect. 
Though the perturbation scales as b, the time for which 
it is important scales as b~ 2 . So Fig. 2c shows that 
b = 0.4, which saturates at ~ 10 kicks, gives < p >~ 
while b = 0.02, which saturates at ~ 5000 kicks, gives a 
very large asymmetry. 

For the standard map/QKR there is a well-known re- 
lation between the quantum localization length and the 
classical diffusion constant: L ~ where the constant 
a was found to be 1/2 Q. The N = 3, b = 0.1 cdw-QKR 
takes a modified proportionality constant, ie L ~ *Mp. 
We have found this by fitting DL forms to a set of quan- 
tum distributions N(ja) corresponding to different D and 
ft. For large p, the distribution is well approximated by 
the DL exponentially localized profile. 

In Fig. 3a we plot a set of calculated L (which range 
from L ~ 10 — 80) against D for ft = 1/2 together with 
the line corresponding to L = ^4^- The agreement is 
excellent. From L 2 ~ Dt* we obtain t* ~ 12D/H 2 . 

In Fig. 3b, for the quantum distributions in Fig 3a, we 
have also plotted the net momenta as a functions of D, 
together with their classical equivalents, obtained from 
an ensemble of 500,000 classical particles. We see that 
the classical < p > fall monotonically with D, apart from 
fine structure which we attribute to cantori. The quan- 
tum results however, for low D, are much smaller than 
the classical values but increase in magnitude until there 
is a 'cross-over' point at D ~ 3, after which they are 
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FIG. 3: (a) Relation between classical diffusion rate D and 
the quantum localization length L for the cdw-QKR. The solid 
line corresponds to L = 3 'S . (b) Average momentum plotted 
against D for a quantum wavepacket (crosses) and a classi- 
cal 'gaussian wavepacket' (circles). The graph illustrates the 
fact that if the quantum break-time is too short (low D) the 
quantum momentum is small, but catches up with the classi- 
cal momentum at t* ~ t T - 



much closer to the classical values. 



We do not expect perfect agreement with the classical 
results for ft = 1/2; a cleaner comparison might be ob- 
tained for smaller ft, but this might be harder to achieve 
in an experiment. We estimate the quantum break-time 
at the cross-over t* ~ 12D/K 2 ~ 150 kicks. The ratchet 
time t r ~ 5/(Db 2 ) ~ 160 kicks. Such good agreement is 
somewhat fortuitous, since there are larger uncertainties 
in the timescales. Nevertheless it does provide us with a 
useful guide for the best parameters for an experiment. 

So one of our key results is that the requirement t* ~ t r 
implies that we need Db/h ~ 1. This is a necessary 
but not sufficient condition to ensure the best cold-atom 
ratchet. Clearly the L values should be experimentally 
plausible L ~ 10 — 100, so this places a constraint on 
D/h. Equally, for a good classical ratchet, b should be 
small ie b < 0.2. 

Although the values of D considered may be small com- 
pared with those of the standard map, we find the chirped 
dw-QKR is much more chaotic at a given value of D. 
There are no remaining stable islands visible on an SOS 
at the parameters corresponding to Figs 1-3. 

In conclusion, we have found a mechanism for a 
generic, completely chaotic Hamiltonian ratchet. But 
much remains to be done, such as considering the ef- 
fect of the partial loss of coherence in the experiment or 
the effect of the finite duration of the experimental kicks. 
Though transport appears as a general sign-dependent 
correction to the classical diffusion rate , the role of spe- 
cific dynamical effects such as desymmetrization of Levy 
flights pi must be investigated. 

We thank Prof. S Fishman for helpful advice. N.H 
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